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ABSTRACT 


The gamma-ray line produced at 0.51-MeV in solar flares is the 
result either of free annihilation of positrons with electrons 
or of the decay of positronium by 2-photon emission. Positron 
annihilation from the bound state of positronium may also pro- 
ceed by 3-photon emission, resulting in a continuum with energies 
up to 0.51-MeV. Accurate calculations of the rates of free 
annihilation and positronium formation in a solar-flare plasma 
are presented. Estimates of the positronium-formation rates 
by charge exchange and the rates of dissociation and quenching 
are also considered. The temperature and density dependence 
of the ratio of 3-photon to 2-photon emission is obtained* For 
temperatures less than 10 K, this ratio is found to depend pri- 
marily on the density of the annihilation region. The observa- 
bility of the 3-photon emission is hindered by the flare-produced 
continuum radiation and other anticipated line emission but aided 
by the inherent delay in the production and slowing-down time of 
the positrons. Asymmetrically-broadened 0.51-MeV line emission 
could be detected at times late in solar gamma-ray events when 
the continuum and prompt line emission have essentially disap- 
peared. When the ratio of free electrons to neutral atoms in the 
plasma is approximately unity or greater, the Doppler width of 
the 0.51-MeV line is a function of the temperature of the annihi- 
lation region. For the small ion densities characteristic of the 

photosphere, the width is predominantly a function of the density. 

0 

For temperatures greater than 10 K, the rate of positron conver- 
sion into -rays is the free annihilation rate. At lower tem- 
peratures this conversion rate is determined by the rates of 
thermalization and positronium formation, which is faster than 
free annihilation by at least two order of magnitude. 
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1 . INTRODUCTION 

Gamma-ray line emission at an energy of approximately 0.61 
MeV was observed by Chupp et al. (1973,75) from the 1972 August 4 
and August 7 solar flares. This line is believed to be due to 
the annihilation of positrons which result from the decay of 
mesons and radioactive nuclei produced in nuclear reactions of 
flare accelerated particles with constituents of the solar atmos- 
phere (Lingenfelter and Ramaty 1967; Ramaty and Lingenfelter 1973; 
Ramaty, Kozlovsky and Lingenfelter 1975) . 

The formation of the 0.51-MeV line depends on the sources of 
the positrons, on the propagation of the positrons in the solar 
atmosphere, on the density and temperature of the ambient medium 
in which the positrons slow down and annihilate, and on the mode 
of positron annihilation since the positrons may annihilate freely 
or from a ground state of positronium. 

Four distinct observable parameters related to the annihila- 
tion of positrons may yield information about the ambient medium 
or the source of positrons. These are the width of the 0.51-MeV 
line, the strength of this line relative to the intensity of 
OLher gamma lines, the strength of the three gamma continuum below 
0.51-MeV which comes from triplet positronium decay, and the 
time dependence of the 0.51-MeV line. The time dependence has 
already been analyzed by Wang and Ramaty (1975), and the number 
of positron emitters relative to the number of neutrons has been 
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related to the spectrum of primary particles (Ramaty, Kozlovsky, 
and Lingenfelter , 1975). In the present paper we investigate in 
detail the slowing down and annihilation of positrons and the 
formation of positronium in a solar flare plasma. We then con- 
sider how the width of the 0.51-MeV line and its strength rela- 
tive to the 37 ^ continuum from positronium decay depend on the 
temperature and density of the medium in which the positron 
comes to rest. 

Positrons in a solar flare are created with energies of the 
order of 1 MeV or higher (e.g. Ramaty, Kozlovsky and Lingenfelter 
1975 ) . They subsequently lose energy at least as fast as the 
standard ionization energy loss of charged particles in matter. 
Approximate calculations indicate that only 10% or fewer of the 
positrons undergo annihilation in flight before slowing to thermal 
velocities, if they are in a plasma, or before slowing to veloci- 
ties at which positronium formation can occur in the neutral gas. 

In general, the Doppler shift of the gammas from the free annihila- 
tion of high energy positrons will be so great that the gammas will 
not be counted as part of the 0.51-MeV line. 

Positronium, denoted by the symbol Ps , is the bound state of 
a positron and an electron. In a fully ionized plasma it can be 
formed by the two body radiative recombination reaction with free 
electrons, or by three body processes only important at very 
high densities. 
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The positronium atoms are formed in triplet or singlet 
spin states according to the statistical ratio of 3:1. If 


undisturbed for times long compared to their natural decay 


times the positronium decays via the emission of either 3^ 's 




3 3^ with a time constant of ^ 


-7 


-10 


or 2^'s. Ps — - 

sec and ^Ps with a time constant of = 10 sec, for 

triplet and singlet states respectively. The 3/^ decays produce 

a gamma ray continuum with a maximum energy of 0,51-MeV while 

the 2 2^ decays appear as a discrete line. Because the cross 

section for free annihilation of an electron positron pair in a 

triplet state is more than two orders of magnitude smaller than 

for a singlet pair, the observation of a three-gamma continuum 

is proof of positronium formation and decay. 

Positronium once formed in a plasma can be dissociated by 

collisions with free protons or electrons whenever the energy 

available in the center-of-mass frame is above 6.8 eV. Since 

_ 2 

the atomic cross sections are typically of the order of'T^a^ , 

the critical density for dissociation of the triplet state is 
14 -3 

about 10 cm , while the singlet positronium begins to be broken 
17 -3 

up at N ^ 10 cm . Another process which depletes positronium 

in the triplet state is the flipping of the spin from triplet to 

singlet through elastic or inelastic collisions with free electrons. 

The density at which this quenching occurs is also of the order of 
-3 

10 cm . 
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Since temperatures high enough to maintain a fully ionized 

plasma do not ordinarily occur in regions of the solar atomsphere 

14-3 

when densities are as high as 10 cm , dissociation and quenching 
processes in a fully ionized plasma are not considered in the 
present work. 

Another set of processes characterize annihilation of 
positrons in a neutral gas. While the energy loss of positrons 
in a plasma is due to many small angle scatterings with electrons 
in the electron gas, in a neutral medium energy loss takes place 
through ionization and excitation of atoms or ions. Once the 
energy of a positron has dropped to about 75 eV, the cross section 
for positronium formation through charge exchange with atoms be- 
comes significant. In fact it dominates the energy loss cross 
sections below 30 eV. As we shall show, in a neutral medium nearly 
all of the positrons form positronium before dropping below 6.8 eV. 
However, if sufficiently high atomic densities obtain, N^IO cm , 
the longer lived triplet positronium will be dissociated by atomic 
collisions or will be quenched through spin flip collisions with 

unpaired atomic electrons. The latter process depletes even those 
3 

Ps atoms which end up with energy below 6.8 eV, the threshold 
energy for dissociation. 

The effect of the dissociation of the triplet state is to in- 
crease the number of annihilations of the positrons either by 

the formation of ^?s subsequent to break-up or by the free annihi- 
lation of positrons that have fallen to energies less than 6.8 eV. 
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The breakup reaction is a very efficient energy loss mechanisin 
since a positron, which goes through the cycle of charge exchange 
followed by breakup, emerges with less than half of its original 

energy. The spin flip reaction also enhances the 2/* annihilations 

31 17-3 

through conversion of Ps to Ps. At densities above 10 cm 

breakup of ^Ps becomes important, as is observed in a typical 

laboratory situation. Since it is extremely unlikely that solar 

positrons can penetrate to regions of such high density, singlet 

breakup is ignored in the present work. 

Another quenching effect, which could deplete the relative 

fraction of ^Ps, is spin flip caused by large magnetic fields. 

In work reported by Wallace (1960) it was shown that magnetic 

fields in excess of 7 kilogauss are needed to reduce the popula- 
3 

tion of Ps states by approximately 30%, at which level the effect 
due to the magnetic field saturates. Such fields far exceed those 
expected in low-density solar regions where other quenching mechan- 
isms are negligible. Since the effects of magnetic quenching in 
all cases are expected to be small, they have not been considered 
in the present work. 

To summarize, in a fully ionized plasma, positrons thermalize 
before either free annihilating or forming positronium. Positron- 
ium is formed by radiative reconibinat ion . In a neutral gas posi- 
tronium formation through charge exchange takes place well before 
thermal izat ion can occur. The solar atmosphere, however, is 
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neither fully ionized nor completely neutral. In the present 
paper, therefore, we investigate the fate of positrons in a 
partially ionized medium. We limit our detailed calculations 
to a hydrogen plasma and we subsequently consider the perturbing 
effects of heavier ions which are found to be negligible. 

In Section 2 accurate cross sections for the free annihila- 
tion, radiative recombination, and charge exchange processes are 
presented. In Section 3 we evaluate the rates of positronium 
fOi..«ation and the resultant energy distributions of the positro- 
nium atoms in a partially ionized medium. The results of this 
section apply to a medium in which the ambient free electrons have 
a Maxwell-Bol tzmann distribution of finite temperature but the 
density of the medium is sufficiently low such that positronium 
atoms decay without further collisions following their formation. 

In Section 4 we additionally consider the breakup of posi- 
tronium as a function of density. The resultant energy distribu- 
tions of the positronium atoms and the relative numbers of triplet 
to singlet pcsjtronium decays are evaluated for temperatures 
characteristic of the solar photosphere. 

In Section 5 we calculate the width of the 0.51-MeV line and 
its strength relative to the continuum as functions of tempera- 

ture and density in the annihilation region. We also consider the 
observability of the 3 7^ continuum in the presence of other solar 
flare continuum radiations and we summarize our results. 
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2. FREE ANNIHILATION, RADIATIVE RECOMBINATION, 
AND CHARGE EXCHANGE 


Free Annihilation 

The annihilation cross section of positrons? with free elec- 
trons, under the assumption that the initial particles are plane 
waves and that the phase shifts introduced by Coulomb interactions 


are negligible, is given by Heitler (1954) . 


where ff mc^ is the total energy of the positron in the electron's 

rest frame, and r = 2.82 x 10 cm is the classical radius of 

o 

the electron. In the nonrelativistic region where most of the 

annihilations take place. Equation (1) reduces to ^pl ane wave, 

2 

non-relativistic) = r c/u, where u is the positron's velocity 

o 

relative to the free electron. The rate of annihilation in this 
case is independent of energy and given by 


where N is the density of electrons, 
e 

Coulomb corrections are important only in the nonrelativistic 
region. When they are taken into account (Landau and Lifshitz 
1958) , the nonrelativistic annihilation cross section can be 


written as 


(T{ VL) 


ZTCoC 


(3) 


= 7C 'Ci 


tA.^ 


/- eyp (-27t<f^) 


where oC is the fine structure constant. Using this cross 
section, t)ie annihilation rate of positrons in thermal equili- 
brium with iu\ isotropic electron plasma of temperature T is 
given by 

U<r(u.) (4) 


where“^j^ and 17^ are the velocity vectors of arbitrary electrons 

3 /2 2 

and positrons, and f ("^ = (m/2i|kT) exp(-mv /2kT) is the 
Maxwel 1-Bolt zmann distribution for particles of mass, m, in 
thermal equilibrium. Equation (4) can be reduced to a single 
integral by transforming the variables of integration to 
- (^) and'TT = • We obtain (Huang 1963) 



i-cjcp ) 


(5) 


The rate coefficients, given by ^ (T)/N , have been evalu- 

£ a e 

ated by numerical integration of Equation (5) and are presented 
in Table 1. As is shown in Figure 1, the effects of Coulomb 
attraction arc most significant at low temperatures. The corrected 
f ree-ann ihi lation rate asymptotically approaches the constant 
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approximate value given by Heitler for temperatures ^ 


TABLE 1. Freo-Annihilation Rate Coefficients 


T in K 

.4 


10 

4 


2 X 10 


6 :( 10 


10 ' 

5 


2 X 10 


6 X 10' 


10 

6 


2 X 10 


4 X 10 


9 X 10 


-14 3 

Rate Coefficient in 10 cm /sec 


14.9 

10.5 

6.1 

4.7 

3.4 

2.1 

1.8 

1.4 
1.2 
1.0 


Positronium Formation by Radiative Recombination 

Positronium formation by radiative recombination is 
essentially the same proce 's as radiative recombination of 
hydrogen except for the difference in the positron and proton 
rest masses. Seaton (1959) has calculated the capture rates 
for hydrogenic ions, based on an as^TTipotot ic expansion de- 
rived by Menzel and Pekeris (1935) and corrected by Burgpss 
(1958) . In the present work, the formulation 


and numerical calculations of Seaton (1959) are used to determine 


the temperature-dependent positron ium-format ion rate. From 
Equations (1) through (7) of Seaton, the total radiative-capture 
rate coefficient is given by 

OC = — J“7 ( m kT J ( ry-f /> 

h 3 i3 m n • 

for hydrogenic ions where the are the coefficients containing 
the contributions from all the angular momentum states for each 
principal quantum number, n. 

According to Seaton (1959), the asympototic expansion enables 
the rates to be calculated with errors not exceeding 2 % for tem- 

4 

peratures of order 10 K, but possibly as great as 20% for tempera- 
tures of order lO^K or greater. It shoxild be noted that the 
widely-used approximation derived by Kaplan and Pikel'ner (1970) 
is obtained from a similar expression for the radiative recombina- 
tion rate which is in agreement with that from Seaton (1959). But 
the expression for the cross-section given by Kaplan and Pikel’ner 
(1970) in their Equation (1.21) appears to be in error, that is, 
too low by a factor of two. Since the expression presented ly 
Kaplan and Pikel'ner in Equation (1.22) is in agreement with 
Seaton, it may be assumed that the discrepancy is only a typo- 
graphical error. 

In order to calculate the radiative capture rate for posd- 
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tronium from Equation (6), the appropriate reduced mass, m/2, 

must be substituted for each m. This is equivalent to putting 
3/2 

a factor of 2 in the coefficient before the sum, and selecting 

the C (^) from Seaton's tabulation for which the arguments are 
n 

equal. That is 

^ \ . <■/) 

T 

for positror:. 'n thermal equilibrium. 

The results of this calculation are presented in Table 2, 
and are shown for compa,-ison with the free ...inihilat ion races in 
Figure 1, Also for comparison the radiative capture rate calcu- 
lat 3 d by Nieminen (1967) is shown as a dashed line. Nieminen's 
calculation of the radiative capture rate is performed as a func- 
tion of velocity which he t)ien relates to an oquivalcr*- tempera- 
ture. When Nieminen's values are folded into the Maxwell -Boltzmann 
velocity distribution, the resulting rate coefficient is in complete 
agreement with the temperature-dependent, radiative-capture rate 
presented here. Frcm Figure 1 it can be seen that including the 
effects due to the Coulomb interaction and the Maxwell-Boltzmann 
di8t**ibut ion yields a small change in tiie temperature at which 
the rates are equal, but yields a large increase in the calculated 
rates. 
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TABLE 2 , Radiative Recombination Rate Coefficients 

* 3 

T in K C Rate Coefficient in cm /sec 

n 


5 X 10^ 

1.998 

1.65 

X 

10-1^ 

2.5 X lo'^ 

1.324 

4.89 

X 

10-“ 

5 X lo"^ 

1.055 

2.76 

X 

10-“ 

2.5 X 10^ 

0.529 

6.10 

X 

10-“ 

5 X 10^ 

0.364 

3.00 

X 

10-“ 

2.5 X 10^ 

0.127 

4.68 

X 

10-“ 


*In Seaton's notation, C = n S ( A )• 

n — n 
n 

From Figure 1 we conclude that in a fully ionized hydrogen 
plasma f ee annihilation and radiative vtrcombinat ion are of 
comparable size at T 8 x lO^K. At higher temperature's Tree 
annihilation dominates and at lower temperatures radiative com- 
bination becomes the dominant process. 




Positron ium Formation by Charge Exchange 

The alternate reaction producing positronium is the charge 
exchange of an electron by a positron with atoms or ions in the 
solar atmosphere. Charge exchange with heavy ions is always 
suppressed by their low number den':i’ly relative to hydrogen. At 
high temperai.ures, T > lO^K, this reaction is additionally sup- 
pressed by the high ionization potentials of the residual ions. 
Estimates indicate that this process cannot compete with the free 
annihilation ri-:es presented previously. At lower temperatures, 

T ^ 10 K, the only effect of charge exchange with heavy ions is to 
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increase slightly the rate of positronium formation. Therefore, 

for simplicity, we consider charge exchange with hydrogen only. 

Published cross-sections for charge exchange with neutral 

hydrogen vary greatly, ranging from peak values of 0,065 (Pels 

2 

and Middleman 1967) to 3,6 (Cheshire 1964) in units of 71" a , 

o 

where a is the Bohr radius, A recent calculation by Chan and 
o 

Fraser (1973) obtains an S-wave contribution with a peak value 
-3 2 

of 3 X 10 7T a , for an energy slightly above th.e threshold 
o 

energy of 6.8 This value, obtained with a difficult but 

presumably accurate variational technique is significantly 
lower than that expected from partial wave calculations employ- 
ing the Born approxim.ation. As a best estimate of the charge 
exchange cross section of positrons with neutral hydrogen, we 
have used a calculation by Drachman et al. (to be published) 
which uses the Born approximation for all partial waves besides 
the S and P. The S-wave result of Chan and Fraser is used and 
the P-wave is modified by scaling the Born partial wave to 
agree with that calculated by Chan and McEachran (to be pub- 
lished) at one energy just above threshold. The resultant 
cross section rises more slowly than the distorted-wave calcula- 
tion according to Massey and Mohr (1954) but peaks at approxi- 
— 2 

mately 3.3 at about 14 ''' positron energy. The values of 

the cross section for charge exohi^ge with hydrogen as calcu- 
lated by Drachman et al. and used in the present work are given 
in Table 3. 
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TABLE 3. Charge Exchange Cross Sections on Atomic Hydrogen 
Positron Energy Cross Section 


(eV) 

<71 a 

7.5 

0.50 

10.0 

2.4 

15.0 

3.3 

25.0 

2.0 

50.0 

0.41 

75.0 

0.12 


Detailed calculations of rate coefficients based on the 
cross sections given in Table 3 are complicated by the possi- 
bility that positrons may form positronium before they therma- 
lize, or that once thermalized, their energies may be below 
the positronium-formation threshold. These questions are 
treated in the next section. Nevertheless an estimate of the 
significance of the charge-exchange process may be obtained 

by noting that at 10 eV, for example, the charge exchange rate 

-13 3 -1 

coefficient is approximately 8 x 10 cm sec for a relative 

5 -5 

neutral hydrogen abundance at 10 K of 2 x 10 (Gabriel 1971) . 
This rate coefficient is larger by about an order ot magnitude 
than the radiative-recombination rate coefficient shown in 


Figure 1 at this temperature. 
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3. POSITRONIUM FORMATION AND ANNIHILATION IN A LOW DENSITY PLASMA 

The calculations presented in the previous section demonstrate 

that the fate of positrons in a fully ionized plasma at densities 
14 -3 

less than 10 cm is completely determined by the ratio of the 
thermally averaged radiative recombination rate to the free anni- 
hilation rate. However, because positrons result from hadronic 
reactions with relatively long interaction lengths, a more real- 
istic assumption may be that the annihilations occur in the chromosphere 

or photosphere. Such considerations necessitate the investigation 

4 

of positron interactions in media with temperatures as low as 10 K 
or somewhat less and with a wide range of densities and ion frac- 
tions. In this section we re-examine the assumptions implicit 
in the application of thermally averaged rates to determine the 
fate of positrons in solar flares. Calculations are presented 
showing to what extent positrons thermalize before annihilating 
or forming positronium. For those positrons whose energy falls 
below 6.8 eV, the threshold for charge exchange with a hydrogen 
atom, the competing processes of free annihilation and thermal- 
izing up (populating the high energy tail of a thermal distribu- 
tion) are evaluated. To facilitate these calculations the effects 
of dissociation and quenching are deferred to the Monte Carlo 
treatment described in Section 4. 

Energetic charged particles, which enter a plasma of lower 
chara'^^-erist ic thermal energy, slow down due to long range Coulomb 
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interactions. One method of treating the energy loss is by use 
of a Fokker-Planck equation. The Fokker-Planck equation neglects 
the generation of plasma waves but includes the Coulomb interac- 
tion. Neglect of plasma oscillations in a thermal plasma under- 
estimates the energy loss rate by about 20%. See, for example, 

Sigmar and Joyce (1971) . In order to calculate the relative 
rates of positronium formation, terms which represent charge 
exchange with neutral hydrogen and radiative capture must be 
added to the Fokker-Planck terms. The resulting equation is 

‘ fPfMt)} (Ne^r 

( 8 ) 


where f (\T, t) is the distribution function of the particles 
entering the plasma; fP 1 1 is the Fokker-Planck operator 
(Rosenbluth et al. 1957, Montgomery and Tidman 1964); N^ is 
the density of neutral atoms; is the plasma density; 

is the radiative capture cross section (Nieminen 1967); 

•'rr 

and CT_ is the charge exchange cross section from Table 3. 
ce 

Details of this calculation have been presented by Joyce et al 

(University of Iowa Research Report 76-2, unpublished). 

In this treatment, the energy loss due to collisions with 

the neutral gas is neglected which is a reasonable approximation 

if N /N 0.1 . The two cross sections for positronium forma- 

e n 

tion dominate the positron loss rates for the range of temperatures 
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and ion densities considered. We have considered plasmas with 

the following properties; the temperature of the ambient plasma 

4 

and neutral gas is 1.16 x 10 K corresponding to a kinetic energy 

of 1 eV. The number density of the neutral gas plus the plasma 
12 -3 14 -3 

is 10 cm or 10 cm . The ratios of plasma to neutral gas 
are 1 or 0.1. The initial distribution function of the positrons, 
f(^o), is a delta function corresponding to an energy of 50 eV. 

The results of the integrations of Equation (8) are summarized in 
Figures 2 through 7. 

Figure 2 shows the evolution of the positron velocity distri- 
bution, N(v,t) = 4iTv f(^t) for equal ion and neutral densities. 

The total number of positrons has been decreased by about 28% by 
the time the distribution is almost thermal (tfe 3 x 10 ^ sec). 

Since the charge exchange cross section drops to zero at a thres- 
hold energy of 6.8 eV, positronium formation continues quite 
slowly after the bulk of the particles have energies characteristic 
of the plasma thermal energy. 

Figure 3 shows the positron evolution for = 0.1 and 

14 

N + N =10 . In this case 94% of the positrons form positronium 

e n r IT 

by charge exchange before the distribution thermalizes. 

Figure 4 shows the fraction of positronium having formed as 
a function of time for two ratios of plasma to neutral density. 

There are two regions of development. First, the positronium 
fraction rises rapidly until the remaining positrons have thermalized. 


This is followed by a region of slow increase which is approximately 
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1 inear in time. By selectively excluding some of the processes 
in the equation, we have determined that the slow rise is due 
to charge exchange of positrons in the thermal tail of the dis- 
tribution which extends to the energy region above the charge 
exchange threshold. As positrons in the thermal tall are lost 
through positronium formation, the tail is repopulated by posi- 
trons which gain energy. 

For N /N =0.1, most of positronium is formed in the first 
e n 

stage so that the second stage is of no great importance. The 

duration of the first stage is approximately 10 ^ sec, consistent 

with the charge exchange cross section given in Table 3, 

For N /N = 1.0, a large fraction of the positronium is 
e n 

formed in the second stage. Since the number of positronium 

atoms exhibits linear grovv _h, we have extrapolated the curve 

for N /N = 1.0 to determine the time for which all the positrons 
e n 

can be converted to positronium. The results of this extrapola- 
tion are presented in Table 4. 


TABLE 4. Time to Convert all Positrons of Initial Energy 50 eV 

to Positronium in a 1-oV Plasma with N = N , 

e n 


N + N 
e n 

, -3 

in cm 


Positronium Formation 
Time in sec 


CM 

O 

2.3 

X 

10-2 

io“ 

2.4 

X 

10-2 


to 

00 

X 

10-^ 
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The characteristic times to annihilate by free annihilation 

or to form positronium by radiative recombination in a medium 

of any specific density and temperature can be calculated from 

the corresponding rate coefficients presented in Tables 1 and 2. 

4 

A comparison of these times for a temperature' of 1.16 x 10 K 
and the densities considered in Table 4 shows that the charge 
exchange process, including the time to thermal ize up, dominates 
over the competing processes of free annihilation and radiative 
recombination. 

Figures 3 and 6 show the distribution of energies at which 
positronium has formed as a function of time. At large times a 
sharp peak begins to develop at low energies. This is due to 
positrons being thermal ized upward until they are above the thres- 
hold energy for the charge exchange cross section. This effect 
can be seen more clearly if the plasma temperature is slightly 
larger. Figure 7 was calculated for a plasma with kT = 4 eV, 

Since more of the thermal tail overlaps the charge exchange thres- 
hold, the development of the peak occurs on a shorter time scale 
than in the 1 eV plasma. In all of these calculations, the energy 
loss due to elastic and inelastic collisions with neutrals has 

been neglected. For the case N = N this is a good approximation, 

e n 

but for smaller ionization fractions, the neutral particle interac- 
tions can cause a significant amount of energy loss. The effect of 


including energy loss on neutrals and the generation of plasma 
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waves, mentioned previously would result in a faster energy loss 
so that the mean energy of positronium at formation would be 
less than predicted from these calculations. 

The results of the calculations presented in this section 

are applicable to positron interactions in media for which the 

sum of the electron density and the neutral ion density, N = N + N , 

e n 

14 -3 

is bounded by N;^ 10 cm . For the chromosphere, this density 
is a reasonable upper limit and justifies the neglect of disso- 
ciation of positronium by collisions with electrons, protons, 
or hydrogen atoms. On the other hand, the observed rise time of 
the 0.51-MeV line implies densities, N ^ 10 cm in the annihi- 
lation region (Wang and Ramaty 1975) . In the chromosphere, den- 
sities in the range 10^^ to lO"*"^ cm ^ have been inferred from 
analyses of white-light flares (Hudson 1972, Machado and Rust 1974) . 
However, because the rise-time observations set only lower limits 
on the density, annihilation in denser regions such as the photo- 
sphere is also possible. In determining the parameters associated 
with positronium formation in media for which N ^ 10 cm , the 
effects of dissociation and quenching may not be neglected. In 
the next section these effects as well as approximations for the 


energy loss processes with neutral atoms and plasma are considered. 
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4. POSITRONIUM FORMATION, ANNIHILATION, AND DISSOCIATION 
IN A LOW TEMPERATURE PLASMA 

In this section we present a simplified Monte-Carlo calcu- 
lation of the fate of a positron annihilating in a predominantly 

neutral gas of temperatures kT 1 eV and density in the range 

12 17 -3 12 -3 

10 N ^ 10 cm . The density of 10 cm is inferred 

from the rise time of the 0.51 MeV line as discussed in the 

17 -3 

previous section. Densities up to 10 cm correspond to photo- 
spheric annihilation regions. Indeed, in the exposition of their 
model of white light flares, Najita and Orrall (1970) show that 
energetic ions between 10 MeV and 1 GeV release most of their 
energy to the ambient gas in the photosphere. As far as posi- 
tron reactions are concerned, that region is more nearly a 
neutral atomic gas than a plasma. 

The principal mechanism for energy loss by positrons inter- 
acting with neutral constituents of the ambient medium is the 
ionization or excitation of hydrogen atoms. Omidvar (1965) has 
shown that in the ionization process, the electron carries off 
a kinetic energy equal, on average, to h its binding energy of 
13.6 eV. This means that the positron loses about 17 eV in each 
ionizing collision. Since the lS-^-2P excitation costs 10.2 eV, 
we have taken the average loss to be about 13.6 eV and the energy 
loss cross section to be the sum of the separate cross sections 
for ionization and excitation of hydrogen. 


In addition to the neutral atoms there are also free elec- 


trons and ions in the solar atmosphere. As we saw in Section 2 , 
the scattering of positrons from free electrons is a very effi- 
cient energy loss process. To compare it with the atomic process 
we can use the expression listed by Green and Wyatt (1965) . For 
the effective cross section for the energy loss E by positrons 
of energy E in a plasma, we take 

cr = sn:e'^_ ^ ( Xo /h) ■ 

While the effective cross section is written as if it were for 
a single scattering event, it is actually built up from many 
small angle scatterings. ^ 0 = the Debye length, the range of 
Coulomb force being made finite by shielding, and b is the dis- 
tance of closest approach. One finds that for 50 eV positrons 
with ^ E = 13.6 eV, for example, the effective cross section is 
two order of magnitude larger than the cross sections for ioniza- 

_3 

tion and excitation of hydrogen. Thus for N^/N ^ 10 the free 
electrons cannot be ignored. Since the standard solar atmosphere 
with lO^^cm ^ N lO^^cm ^ comprises relative electron densi- 
ties of this magnitude, we took into account loss against the 
free electrons by treating the cross section of Equation (9) in 
the same way as the ionization and excitation cross sections. 

Our simplfied Monte Carlo calculation has been carried 
out by assigning each process an average energy loss so that 
only a small set of positron or positronium kinetic energies 
are associated with each initial positron energy. Then the 
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probability for either triplet or singlet positronium decay 
can be written as a finite series of terms, each term being 
equal to the probability of occurrence for the sequence of re- 
actions it represents. The four processes are: energy loss 

to neutral hydrogen atoms or free electrons, /k E • 13.6 eV; 
positronium formation through charge exchange, A E ■ 6.8 eV; 
breakup of the positronium through collisions with hydrogen 
atoms, A E * (E (Ps) + 6.8)/2 eV; and spin flip of triplet to 
singlet, A E « 0. 

From the Monte Carlo calculations we obtain the probabili- 
ties that a positron decays via triplet positronium, P(^Ps— ^3j^ ) 
or via singlet positronium, P(^Ps-^2^). A positron that falls 
below 6.8 eV cannot form positronium so that only free annihila- 
tion is possible unless the ion density is sufficient for it to 
thermalize up (see Section 3). The probab.Hity of free annihila- 
tion is calculated from 

P(f3) ~ 1 - rCPs^sy) . 

The rates for each process, except the triplet decay rate, depend 
on density and energy through the relation 

\ - Alu cr 

where N is the number de'.'^ity of target particles, u is the 
relative velocity of the reactants, and U* is the total cross 


( 10 ) 


( 11 ) 
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section for that process. As previously noted, for densities 
characteristic of the solar atmosphere, only the triplet posi- 
tron ium will be broken up or quenched since the singlet decay 
rate is large compared to the rates of the alternative atomic 
processes. 

The cross sections used in the calculation are taken either 
from electron scattering experiments on hydrogen atoms or from 
theoretical calculations. The cross section for the ionization 
of hydrogen by positrons is taken to be equal to the electron 
ionization cross sections reported in Golden and McGuire (1974) . 
The cross section for the 1S-»>2P excitation of hydrogen by posi- 
trons was taken to be the same as t' experimental electron cross 
section to be found in Geltman (1969) . The breakup cross section 
was taken from Massey and Mohr (1954) while the spin flip cross 
section is a new result of Hara and Fraser (1975) . The charge 
exchange cross section is from Table 3. The temperature is 
assumed to be so low that the thermal motions in the ambient 
medium are negligible relative to the kinetic energies of the 
incident p:j.,itrons. The relative ion densities are taken from 
the standard solar atmosphere according to Allen (1973) for each 
total density considered in these calculations. 

Our results are displayed in Figure 8. At low densities, 

the positrons annihilate through pos/.^ronium decay with a 3:1 

triplet to singlet ratio. As the density increases through 
15 -3 

10 cm the triplet positronium is either broken up or quenched 
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through rpin flip. In either case mostly singlet positronium 
is produced although the number of free positrons ‘/ith energies 
less than 6.8 eV increases. The probability of free annihila- 
tion has probably been ovi*x estimated since we have ignored 
those positrons which replenish the high-energy thermal tail 
above 6.8 eV. However, this error does not affect the estimate 
of the number of 2^ photons per positron, * since 

12 -3 17 -3 

^2"^ varies between 0.5 for N = 10 cm and 2.0 for N = 10 cm . 

At the lower densities the ^Ps is formed, on average, at 
higher kinetic energies since it does not arise from the breakup 
of ^Ps and the subsequent ^Ps formation. Since the decay rate 
is too rapid for thermal izat ion of the potitronium, our Monte Carlo 
calculation also enables us to find the average ^Ps kinetic energy 
before decay. The average energy is sl'own in Figure 9, from 
which it is clear that the mean energy decreases with increasing 
density. 

If we assume that positrons which fall below 6.8 eV gradually 

thermalize, at T ^ 1 eV, most will free annihilate with the 

+ — 

center of mass motion of the e e pair being determined by the 
Fermi motion of the IS hydrogenic electron. We can assign an 
approximate kinetic energy of 6.8 eV to this motion. The dotted 
curve in Figure 9 shows the average kinetic energy for axl 2^ 


sources . 
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5. DISCUSSION 

From the results of the calculations presented in the pre- 
ceding sections, we can describe the fate of a positron in a 
solar flare plasma. 

For temperatures greater than lO^K, positrons annihilate 

-14 -1 

via free annihilation at a rate of about 10 N sec as can 

e 

be seen in Figure 1, For this case the width of the 0.51 MeV 
line is determined by the center of mass motion of the e^ e 
pair which is a function of the temperature of the medium. 

The full-width at half maximum of this line, due to thermal 
broadening, is given by 

^ = 1.1 ^ (kev) , 

where is the temperature in units of lO^K. For T > 10°K, 

A > 11.0 keV. But because solar flare plasmas are not 

expected to be hotter than 3 x lO^K (Neupert, 1968) , 

A < 35 keV. 

For temperatures just below lO^K, radiative recombination 
dominates over free annihilation. The relative importance of 
radiative recombination and charge exchange is determined by 
the residual neutral hydrogen abundance. For the fractional 
neutral densities characteristic of the quiescent solar atmos- 
phere (Gabriel 1971), charge exchange is expected to be the 


(13) 


dominant process. 
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At temperature below a few times 10 K, positronium forma- 
tion through charge exchange is the dominant reaction through 
which positrons annihilate. The relative rates of the competing 
processes of slowing down and of forming positronium depend pri- 
marily on the fractional ion density of the annihilation region. 

When N /N ^0.5, most of the positrons first thermalize 
e n 

and then form positronium by charge exchange. This sequence of 

events could be maintained even if the mean evergy of the ambient 

electrons were less than the threshold for positronium formation 

(6.8 eV) . This result follows from the fact, demonstrated in 

Section 3, that the rate of upward thermalization is larger than 

the rates for both free annihilation or radiative recombination. 

The annihilation rate is determined by the rate of charge exchange 
—8 “1 

(^3 X 10 N^sec , see Table 3) if the mean energy is larger 

thHJ. eV, and by the upward thermalization rate 

( ->10 N sec see Table 4), if this tinergy is less than 
e 

approximately 6.8 eV. Both these rates are much larger than 
the free annihilation and radiative recombination rates shown 
in Figure 1. When the positrons thermalize before forming posi- 
tronium and the mean energy of the ambient electrons is above 
6.8 eV, the width of the 0.51 MeV line is approximately 
^ ^ (T^ - 7) . When the mean energy is less than 

6.8 eV, most of the positrons form positronium at an energy 
just above threshold. The resultant energy of the positronium 
is low, of the order of 1 eV. In this case, ^ 1 keV. 
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When the fraction of ions is low, N /N 0.5, most of 

e n 

the positrons will form positronium before they are slowed 

below akinetic energy of about 15 to 20 eV for which the charge 

exchange cross section is a maximum. The rate of annihilation 

in this case is determined by charge exchange and is approxi- 
“8 “1 A 

mately 3 x 10 sec . Since the decay occurs well 

before the positronium atom can thermal ize, the Doppler broad- 
ening is given by 

= 0.97 (Mean Kinetic Energy in eV) (keV) . (14) 


15 -3 

For high densities, N ^ 10 cm , the mean kinetic energy is 
approximately 13 eV and hence H Ey ^ 3 eV. Also, breakup 
of ^Ps results in no 3 ^ decays. For low densities, N < lO^^cm 


the mean kinetic energy is about 25 eV so that ^ ^ 5 keV. 

(See Figure 9.) 


The values of the measurable quantities 
ratio of the number of 3 3^ to the number of 2^ 


/N^ ^ , the 

decays, and 


^ that may be expected for various regimes of temperatures 

and densities are shown in Table 6. For high temperatures, 

T ^ lO^K, or high densities, N ^ 10^^ cm 

Since the temperature in regions of densities greater than 
14 -3 

10 cm , even during flares, is not expected to be greater 

4 

than a few times 10 K, these two cases can be distinguished by 


the width of the 0.51 MeV line. For high temperatures, 
yt ^ 10 keV; and for high densities, ^ E ^ 


3 keV. 
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For temperatures lower than 10 K and densities less than 
14 -3 

10 cm , all positrons form positronium which annihilates be- 


fore it dissociates. Thus N. 


'3 T ^^2 T 

ionization of the medium and its temperature determine the width 
of the 0.51-MeV line. The largest width, ^ ^ 10 keV is 

obtained for T ^ 10 K, and the smallest width, A St 1 keV 

4 

is obtained for T V x 10 K and > 0*5. However, if 

N /N < 0.5, the width becomes larger, rather than smaller, 
e n 

The annihilation radiation from decay of the triplet state 
of positronium could possibly be observed in measurements of 
the gamma-ray spectrum at energies just below 0.511 MeV with 
detectors having good energy resolution. Such measurements, 
however, are severely complicated by the existence of strong con- 
tinuum emission (Crannell, Ramaty and Werntz 1975). A favorable 
condition for observing positronium annihilation radiation may 
arise at the late stages of solar gamma-ray events when the 
continuum is greatly reduced. This follows from the delayed 

nature of positron annihilation radiation caused by the long 

11 13 

half lives of some of the positron emitters ( C and N) , and 
possibly also by the long slowing down times of relativistic 
positrons from 7T’ ^ decay (Wang and Ramaty 1975). Thus, when 
the numer of accelerated particles in the flare region is 
already diminished and hence no nuclear reactions and brems- 
strahlung are produced, positronium annihilation radiation 
could be more easily observed. 


■4.5. The state of 


Ps DECAY 


THERMALIZATION 


FREE ANNIHILATION 


^3/ / N2y ® 0 

Ey = 1.1 keV 


THERMALIZATION 


CHARGE EXCHANGE 


AEyW(T^-7)^* keV 
Ng/Np > 0.5 




Ne/Np< 0.5 


PARTIAL THERMALIZATION 

I 

CHARGE EXCHANGE 


A E/ 

I 

I 

I 



3ps BREAKUP 







31 


7 . ACKNOWLEDGEMENTS 


The authors are indebted to Dr. R.J. Drachman for stimu- 
lating suggestions. We also wish to thank Drs. R.J. Drachman 
and K.J.H. Phillips for helpful discussions and for suggesting 
references; Drs. R.J. Drachman and K. Omidvar for making the 
results of their calculations available prior to publication; 
Professor S.M. Austin and his colleagues for making their cross 
section measurements available prior to publication (King et al 
1975); and Professor E.L. Chupp and Dr. P.P. Dunphy and for 
making their data available prior to publication. 



32 


REFERENCES 


Allen, C.W, 1973, Astrophys ical Quantities, (Third Edition) Tlie 
Athlone Press, University of London, London. 

Burgess, A. 1958, Monthly Notices of the Royal Astronomical 
Society, 118 , 477. 

Chan, Y. and Fraser, P.A. 1973, J. Phys, B. 2504. 

Chan, Y. and McEachran, R.P. (to be published) . 

Chesire, I.M. 1964, Proc. Phys. Soc. 83 , 227. 

Chupp, E.L., Forrest, D.J., Higbie, P.R., Suri, A.N., Tsai, C., 
and Dunphy, P.P. 1973, Nature, 241 , 333. 

Cnupp, E.L., Forrest, D.J., and Suri, A.N. 1975, Proc. of the 
AU COSPAR Symposium No. 68, Solar Gamma-, X-, and EUV Radiations, 
edited by S. Kane, P. 341. 

Crannell, C.J., Ramaty, R., and Werntz, C. 1975, Proceedings of 
the 14th International Cosmic Ray Conference, Munich, 1656. 

Drachman, R.J., Omidvar, K., and McGuire, J.H. (to be published). 

Fels, M.F. and Middleman, M.H. 1967, Phys. Rev, 163 , 129. 

Gabriel, A.H. 1971, Solar Physics, 21 , 392. 

Geltman, S. 1969, Topics in Atomic Collision Theory , Academic 
Press, New York and London. 

Golden, J.E. and McGuire, J.H. 1974, Phys. Rev. Letts., 32 , 1218. 

Green, A.E.S. and Wyatt, P.J. 1965, Atomic and Space Physics , 
Addison-Wesley Publishing Company, Inc., Reading, Massachusetts. 

Kara, S. and Fraser, P.A. 1975, Journal of Physics B (in press). 

Heitler, W. 1954, The Quantum Theory of Radiation , Oxford Univer- 
sity Press, London. 

Huang, K. 1963, Statistical Mechanics , Wiley and Sons, New York. 
Hudson, H.S. 1972, Solar Physics, 24 , 414. 

Kaplan, S.A. and Pikel'ner, S.B. 1970, The Interstellar Medium , 
Harvard University Press, Cambridge, Mass., USA. 


33 


King, C.H., Rossner, H.H., Austin, S.M. , Chieu, W.S., Mathews, G.J., 
Viola, Jr., V.E., and Clark, R.G. 1975, Phys, Rev, Letts. 988. 

Landau, L.D, and Lifshitz, E.M. 1958, Quantum Mechanics , Pergamon 
Press, Lt d . , London . 

Lingenfelter , R.E. and Ramaty, R. 1967, High Energy Nuclear Reac- 
tions in Astrophysics, edited by B.S.P, Shen W.A. Benjamin, 

New York , 99. 


Machado, M.E, and Rust, D.M. 1974, Solar Physics, 2§.» 499. 

Massey, H.S. and Mohr, C.B.O. 1954, Proc. Phys. Soc. A, 67 , 695. 


Menzel, D.H. and Pekeris, C.L. 1935, Monthly Notices of the 
Royal Astronomical Society, 96 . 77. 

Montgomery, D. and Tidman, D.A. 1964, Plasma Kinetic Theory 
McGraw-Hill, New York, Ch. 2. 

Najita, K. and Orrall, F.Q. 1970, Solar Physics 15 , 176. 
Neupert, W.M. 1968, Ap. J. Letters 153 , L59. 

Nieminen, I, 1967, Arkiv fur Fysik 35 , 1. 

Omidvar, K. 1965, Phys. Rev., 140 , 26. 


Ramaty, R. , Koslovsky, B., and Lingenfelter, R.E. 1975, Space 
Science Reviews, in press. 

Ramaty, R. and Lingenfelter, R.E. 1973, High Energy Phenomena on 
the Sun, edited by R. Ramaty and R.G. Stone, NASA SP-342 (National 
Aeronautics and Space Administration, Washington, D.C.), 301. 


Rosenbluth, N.M., MacDonald, W.M., and Judd, D.L. 1957, Phys. 
Rev. 109, 1. 


Seaton, M.J. 1959, Monthly Notices of the Royal Astronomical 
Society, 119 , 81. 

Sigmar, D.J. and Joyce, G. 1971, Nuclear Fusion, 11 , 447. 

Wang, H.T. and Ramaty, R. 1975, Ap. J., 202 , 532. 

Wallace, P.R. 1960, Solid State Physics , 10 , 1, edited by F. Seitz 
and D. Turnbull, Academic Press, New York. 


34 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


The rates for radiative capture and for free annihi- 
lation of positrons as a function of temperature. 

The solid lines are from the results of the present 
work. The dashed lines indicate previously published 
results in which significant physical considerations 
are neglected. 

Positron distribution as a function of velocity in 

units of v/v^, where v^ is the velocity of an electron 

with kinetic energy kT, at various times for 

14 -3 

N =N =0.5x10 cm and kT = 1 eV. 
e n 

Positron distribution as a function of velocity in 

units of v/v , where v is the velocity of an electron 
e e 

with kinetic energy kT, at various times for 

N =0.1N,N +N = 10^^ cm"^, and kT = 1 eV. 
e n e n 

Fraction of positrons which have formed positronium 

as a function of time for 

14 -3 

N =N =0.5x10 cm and kT = 1 eV; 

e n 

and for N = 0.1 N , 
e n 

N + N = 10^^ cm ^ , and kT = 1 eV. 
e n 

Fraction of positrons which have formed positronium 

as a function of energy at various times for 

14 -3 

N =N =0.5x10 cm and kT = 1 eV. 

e n 

Fraction of positrons which have formed positronium 

as a function of energy at various times for 

N =0.1N ,N +N = 10^^ cm’^, 
e n e n 


and kT = 1 eV. 
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Figure 7. 


Figure 8. 


Figure 9, 


Fraction of positrons which have formed positronium 

as a function of energy at various times for 

14 -3 

N *N *0.5x10 cm and kT = 4 eV. 
e n 

Probabilities of 3 and 2/* annihilation of posi- 
tronium as a function of density in the standard solar 
atmosphere. The fraction of free annihilation may be 
overestimated since thermal izat ion upward is not con- 
sidered. 

Average kinetic energy of the center-of-mass of the 
two photons resulting from singlet positronium at 
decay as a function of the density of the ambient 
medium. The dashed line includes the two photons from 


free annihilation 
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